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1.  Introduction 


The  phenomenon  of  breaking  waves  in  the  ocean  occurs  whenever  a 
momentarily  high  crest  reaches  an  unstable  condition.  It  occurs  Inter¬ 
mittently,  and  the  frequency  of  occurrence  depends  on  the  severity  of  the 
sea. 

The  significance  of  Information  on  breaking  waves  cannot  be  over¬ 
emphasized.  First  of  all,  breaking  waves  are  always  associated  with  steep 
waves  that  occur  In  a  given  sea;  therefore,  consideration  of  breaking  waves 
cannot  be  Ignored  In  the  statistical  prediction  of  wave  height.  Unfortunately 
currently  available  statistical  prediction  methods  do  not  consider 
the  concept  of  breaking  waves.  Hence,  the  largest  wave  height  expected  to 
occur  In  a  specified  period  of  time  (the  extreme  wave  height)  In  a  given 
sea,  when  predicted  by  using  the  current  prediction  methods.  Is  unreasonably 
high  and  Is  clearly  beyond  the  breaking  wave  limit.  This  situation  can 
be  rectified  by  Introducing  the  concept  of  breaking  waves  In  the  prediction 
of  wave  height. 

Secondly,  breaking  waves  exert  by  far  the  largest  wave-induced  force 
(often  In  the  form  of  Impacts)  on  marine  systems,  which  may  cause  serious 
safety  problems  and  structural  damage  of  the  systems.  Hence,  It  Is  highly 
desirable  to  consider  the  frequency  of  occurrence  and  severity  of  breaking 
waves  In  evaluating  hydrodynamic  forces  for  the  design  of  marine  systems. 

The  area  addressed  by  the  present  study  Is  that  concerned  with 
Information  on  the  frequency  of  occurrence  of  breaking  waves  In  deep  water. 
Breaking  waves  defined  here  are  not  those  whose  tops  are  blowing  off  by 
the  wind  but  those  whose  breaking  Is  associated  with  steepness;  more 


precisely,  breaking  which  results  from  a  steepness  for  which  the  wave 
height  exceeds  14.2  percent  of  the  Stokes  limiting  wave  length  (Mlchell  1893). 
These  wave  phenomena  are  extremely  complicated  and  the  times  of  their  occurrences 
are  usually  unpredictable.  Hence,  they  may  be  most  usefully  evaluated  through 
the  probabilistic  approach. 

Only  a  few  studies  have  been  carried  out  on  the  frequency  of  occurrence 
of  breaking  waves  In  deep  water.  Nath  and  Ramsey  (1976)  assumed  that  the 
wave  height,  H,  and  Its  period  squared,  T^,  were  statistically  Independent 
and  that  both  follow  the  Rayleigh  probability  law.  On  the  other  hand, 

Houmb  and  Overvlk  (1976)  evaluated  the  probability  of  breaking  using  the 
joint  probability  function  of  wave  height  and  period  applicable  for  a 
narrow-band  random  process. 

Although  the  assumption  of  the  narrow-band  random  process  was  used 
In  both  of  the  above  studies,  the  wave  records  which  Include  very  steep 
waves  whose  breaking  Is  Imminent  Indicate  that  the  assumption  may  be  highly 
unrealistic.  Hence,  In  the  present  study,  a  method  to  evaluate  the  frequency 
of  breaking  waves  Is  developed  based  on  the  joint  probability  distribution 
of  wave  height  and  period  for  a  non-narrow-band  random  process.  The  method 
of  approach  Is  outlined  In  the  following: 

In  general,  for  waves  with  a  non-narrow-band  spectrum.  It  may  be  well 
to  consider  breaking  which  occurs  under  two  different  situations;  one  which 
takes  place  along  an  excursion  above  the  zero-line,  the  other  which  occurs 
along  an  excursion  crossing  the  zero-line.  This  situation  is  shown  in  an 
explanatory  sketch  In  Figure  1.  In  the  figure,  A  and  B  represent  the  positive 
maxima  and  minima,  respectively,  of  a  non-narrow-band  random  process,  while 
C  and  D  represent  the  negative  minima  and  maxima,  respectively.  Note  that 
possible  occurrences  of  breaking  are  expected  on  the  excursion  CA  as  well  as 


BA.  These  excursions  may  be  called  the  Type  I  excursion.  Type  II  excursion, 
respectively.  Hence,  In  order  to  evaluate  the  probability  of  occurrence 
of  breaking,  it  is  necessary  to  consider  the  joint  probability  distribution 
of  the  excursion  and  its  associated  time  interval  between  two  peaks  as  well 
as  the  frequency  of  occurrence  of  each  type  of  excursion. 

Computations  of  the  probability  of  occurrence  of  breaking  are  made  in 
various  sea  severities  using  a  family  of  wave  spectra  developed  from  an 
analysis  of  wave  data  measured  in  the  North  Atlantic  Ocean.  From  the  results 
of  the  computations,  the  limiting  sea  severity  is  obtained  below  which  no 
wave  breaking  is  expected  to  occur.  It  is  also  found  that  the  probability 
of  occurrence  of  breaking  increase  significantly  with  increase  in  the  4th 
moment  of  the  wave  spectrum  irrespective  of  sea  severity. 

2.  Wave  Breaking  Criterion 

Prior  to  developing  a  formulation  for  evaluating  the  frequency  of  occur* 
rence  of  wave  breaking  in  irregular  seas,  it  may  be  well  to  discuss  the 
criterion  of  the  wave  breaking  phenomenon  which  will  be  used  in  the  prediction. 

The  criterion  for  breaking  of  regular  waves  in  deep  water  has  been 
addressed  from  various  viewpoints.  These  include  Stokes  (1880),  Mitchell 
(1893),  Longuett-Higgins  (1963,  1969,  1974,  1976),  Cokelet  (1977),  etc.  The 
most  commonly  known  breaking  criterion  is  that  the  wave  height  should  exceed 
14.2  percent  of  the  Stokes  limiting  wave  length  which  is  20  percent  greater 
than  that  of  ordinary  sinusoidal  waves  of  the  same  frequency.  That  is,  the 
breaking  criterion  is  given  by, 

0.142  (1) 
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where,  H  =  wave  height 

L^=  Stokes  limiting  wave  length  =  1.2g/(2iTf^) 
f  »  frequency  in  Hz 

Hence,  in  terms  of  wave  period,  wave  breaking  occurs  when  the  wave  height- 
period-relationship  becomes, 

H>  0.027  gT^  (2) 

(1968) 

Dean  gives  a  value  of  0.033  for  the  magnitude  of  the  const  in  Equation 

A 

(2)  based  on  his  computational  results  using  the  stream  functio  resentatlon 
of  nonlinear  waves.  On  the  other  hand,  the  results  of  laboratory  tests 
carried  out  by  Van  Dom  &  Pazan  (1975)  show  that  all  steep  regular  waves 
generated  in  deep  water  break  at  some  lower  wave  height  than  that  given  in 
Equation  (2). 

Since  there  are  no  specific  criteria  for  breaking  of  Irregular  waves, 
the  condition  given  in  Equation  (2)  was  first  used  in  this  study  for  evaluating 
the  frequency  of  wave  breaking  occurrences.  However,  the  results  of  tests 
on  the  frequency  of  breaking  of  Irregular  waves  generated  in  the  tank  have 
shown  that  the  observed  number  of  breakings  during  a  specified  time  period 
is  several  times  greater  than  that  theoretically  computed  based  on  the 
criterion  given  in  Equation  (2).  In  order  to  clarify  this  difference, 
an  experimental  study  was  carried  out  on  the  wave  height-period-relationship 
when  irregular  waves  break.  It  was  found  that  the  breaking  criterion  for 
irregular  waves  is  substantially  different  from  that  known  for  regular  waves. 
The  details  of  the  tests  on  the  breaking  criterion  applicable  for  irregular 
waves  are  discussed  below. 

Several  series  of  irregular  waves  with  different  wave  spectra  were 
generated  in  the  40-m  long  tank  in  the  Coastal  and  Oceanographic  Engineering 


Laboratory.  Waves  were  generated  by  using  pre-programmed  tapes  which  yield 
relatively  severe  waves  leading  to  breaking.  A  capacitance-type  wave-height 
probe  was  installed  at  a  fixed  point  in  the  tank,  and  the  height  and  period 
of  incident  waves  which  were  expected  to  break  right  at  that  location  were 
measured.  Waves  which  were  already  broken  before  they  reached  the  probe 
were  excluded  in  establishing  the  breaking  criterion.  Examples  of  incident 
waves  for  which  breaking  was  imminent  are  shown  in  Figure  2.  Included  also 
in  the  figure  are  the  wave  height,  H.  and  period,  T,  defined  for  establishing 
the  criterion. 

The  wave  height-period-relationship  was  obtained  from  measurements  of 
over  40  incident  waves  about  to  break  and  the  results  are  plotted  in  Figure 
3.  As  can  be  seen  in  the  figure,  the  functional  relationship  between  wave 
height  and  period  at  the  time  of  breaking  of  irregular  waves  is  by  and  large 
different  from  that  known  for  regular  waves.  By  drawing  the  average  line  of 
the  observed  results  plotted  in  Figure  3,  the  following  relationship  is  obtained 
as  the  criterion  for  breaking  of  irregular  waves: 

H  >.  0.020  gX^  (3) 

Let  us  express  the  criterion  given  in  Equation  (3)  in  dimensionless 
form.  As  stated  in  the  previous  section,  two  different  types  of  breakings 
are  considered  in  the  present  study.  For  this,  consider  H  in  Equation  (3) 
as  the  Type  I  as  well  as  Type  II  excursions,  denoted  by  i;»  and  T  as  the  time  interval 
Then,  Q  and  T  are  non-diraensionalized  by  dividing  by  / ra^  (where,  m^  = 
area  under  the  wave  spectrum) ,  and  (where  =  average  time  interval 
between  positive  maxima),  respectively.  This  results  in  the  following 
breaking  criterion  in  dimensionless  form: 


where , 


V  =  dimensionless  excursion  =  C/v^®o 

X  =  dimensionless  time  interval  between  positive  maxima  =  T/T^ 
C  =  excursion 

T  =  time  interval  between  positive  mavima 
Tjjj  =  average  time  interval  between  successive  positive  maxima  A 


e  = 


a  = 


Att 


✓  l-e^ 


1+  /l-i 


band-width  parameter  of  the 

the  j-th  moment  of  the  wave 
0.196  for  m-sec-units 


spectrum  =  /  1  - 


m 


0  4 


spectral  density  function 


3.  Probability  of  Occurrence  of  Breaking  Waves 

As  was  shown  in  the  previous  section,  the  breaking  criterion  consists 
of  two  random  variables,  the  dimensionless  excursion,  v,  and  time  interval, 
X.  Hence,  the  probability  of  occurrence  of  wave  breaking  can  be  evaluated 
from  the  joint  probability  density  function  of  excursion  and  associated 


time  interval. 


The  joint  probability  density  function  of  Type  I  excursion  defined  in 
the  Introduction  (CA  in  Figure  1)  and  the  associated  time  interval,  T,  is 
essentially  the  same  as  that  for  the  positive  maxima.  A,  and  the  associated 
time  interval,  which  was  derived  by  researchers  at  the  Centre  National 
Pour  L'Exploitation  des  Oceans  (CNEXO)  (Arhan  et  al.  1976).  Let  us  assume 
the  wave  profile  x(t)  Illustrated  in  Figure  1  to  be  a  stationary  Gaussian 
process  which  has  an  arbitrary  spectral  bandwidth.  It  is  also  assumed 

f 

that  x(t)  has  zero  mean,  and  the  derivatives  with  respect  to  time,  x(t) 

II 

and  x(t),  both  exist  with  probability  one^  By  considering  twice  the  mag¬ 
nitude  of  the  positive  maxima,  the  joint  probability  density  function  can 
be  written  as  follows : 


f  j.(v  ,  X) 


1 _ (a+jT^hI 


32 /2F 


e(l  -  e^) 


X® 


X  e 


I- 


1  (i  +  /i-gQ 


l-e^ 


where,  v  =  ^  /  v'^  with  ?  = 

magnitude  of  the  positive  maxima 


0  ^  V  <  " 
0  <  X  <  oo 


It  should  be  noted  that  Eq.  (5)  includes  the  positive  maxima  belonging 
to  both  Type  I  and  Type  II  excursions.  Hence,  in  evaluating  the  probability 
of  occurrence  of  wave  breaking  associated  with  Type  I  excursions,  the  maxima 


belonging  to  Type  II  excursions  have  to  be  deleted.  This  can  be  done  by 
taking  Into  consideration  the  frequency  of  occurrence  of  each  excursion. 

That  Isy  the  expected  nimber  of  the  positive  maxima  and  minima,  denoted  by 
and  ,  respectively,  are  given  by  (Cartwright  and  Longuet-Hlgglns  1956) 


The  expected  number  of  the  excursion  BA,  denoted  by  ,  Is  equal 

to  N-.  On  the  other  hand,  the  expected  number  of  the  excursion  CA, 

B 

denoted  by  ,  Is  equal  to  the  difference  of  and  ,  which  In  turn 
Is  equal  to  the  expected  number  of  the  positive  zero-crossings.  Then, 
the  expected  number  of  occurrences  of  Type  I  and  Type  II  excursions, 
respectively,  is  given  by. 


N 


CA 


N 


BA 


(7) 


Hence,  the  frequency  of  occurrence  of  Type  I  and  Type  II  excursions, 
denoted  by  Pj  and  p^j  ,  respectively,  becomes. 
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Thus,  the  conditional  probability  of  breaking  waves  given  that  Type  I 
excursion  has  occurred  can  be  evaluated  from  (5)  and  (8)  taking  into  consider¬ 
ation  of  the  breaking  condition  given  in  (4)  .  That  is. 


Fr{Breaking  wave  ]  Type  I  excursion}  *  P 


00  00 


Then,  the  probability  of  breaking  waves  associated  with  Type  I  excursion 
can  be  obtained  by  multiplying  Eq.(9)  by  the  probability  of  occurrence  of 
Type  I  excursion.  By  taking  account  of  the  fact  that  the  sum  of  the  probabi¬ 
lities  of  Type  I  and  Type  11  breaking  waves  is  equal  to  unity  when  all  waves 
have  broken,  we  have. 


Pr{Breaklng  waves.  Type  I  } 


^  rr  ..... 


A)  dvdA 


For  evaluating  the  breaking  waves  associated  with  Type  II  excursions, 
the  joint  probability  density  function  of  the  excursion  BA  and  the  time 
Interval  is  required.  It  is  given  in  a  dimensionless  form  as  follows: 
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The  derivation  of  the  joint  probability  density  function  is  presented 


in  the  Appendix. 

The  probability  of  breaking  waves  can  be  evaluated  from  Eq. (11)  by 
a  procedure  similar  to  that  given  in  Eq.(lO).  That  is, 

p  rr 

Pr{Breaking  waves.  Type  11}“  - ^ -  I  I  f  (v,X)dvdA  (12) 


The  total  probability  of  breaking  waves  is  the  s\jm  of  the  probabilities 
given  in  (10)  and  (12)  . 

4.  Computation  of  Breaking  Waves  in  Various  Sea  Severities 

It  is  of  Interest  to  examine  the  effect  of  sea  severity  on  the  frequency 
of  occurrence  of  breaking  waves.  For  this,  computations  are  carried  out 
using  a  family  of  wave  spectra  for  a  specified  sea  severity.  The  family 
of  wave  spectra  consists  of  eleven  members  including  the  one  which  is  most 
likely  to  occur  for  a  specified  significant  wave  height,  and  it  represents 
a  variety  of  shapes  of  wave  spectra  observed  in  the  North  Atlantic  Ocean 
(Ochl  and  Hubble,  1976). 

Prior  to  presenting  the  results  of  the  computation,  it  may  be  well  to 
show  an  example  of  the  joint  probability  density  function  of  the  excursion  and 
time  Interval  used  for  evaluating  the  probability  of  breaking  waves. 

Figures  A  and  5  show  the  joint  probability  density  functions  for  Type  I  and 
Type  H  excursions,  respectively,  of  waves  of  significant  height  12.2  m  (40  ft) 
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whose  spectrum  is  shown  In  Figure  6.  Included  also  in  each  figure  is  the  line 
indicating  the  breaking  condition  given  in  Eq.  (4).  If  the  wave  excursion 
exceeds  this  line  for  a  given  tine  Interval  (this  region  is  denoted  by  the 
breaking  zone  in  the  figure),  then  breaking  takes  place. 

As  shown  in  Eq.  (10),  the  probability  of  breaking  waves  associated  with 
Type  I  excursion  is  the  product  of  the  volume  of  the  joint  probability  density 
function  in  the  breaking  zone  indicated  in  Figure  4  and  the  square  of  the 
frequency  of  occurrence  of  Type  I  excursion,  which  is  equal  to  40.66  percent 
for  this  example.  On  the  other  hand,  the  probability  of  breaking  waves 
associated  with  Type  II  excursion  is  the  product  of  the  volume  in  the  breaking 
zone  indicated  in  Figure  5  and  the  frequency  of  occurrence  of  Type  n  excursion, 
which  is  equal  to  10.95  percent  for  this  example.  Thus,  it  appears  that  the 
probability  of  breaking  waves  of  Type  II  excursion  is  approximately  26.9  percent 
of  that  associated  with  Type  I  excursion. 

Figure  7  shows  the  results  of  cosqiutatlon  indicating  the  probabilities 
of  occurrence  of  breaking  in  various  sea  severities.  As  can  be  seen  in  the 
figure,  the  probability  of  breaking  increases  significantly  with  increase 
in  sea  severity,  in  general,  and  no  wave  breaking  is  expected  in  seas  of 
significant  wave  height  less  than  approximately  4  meters  (13.1  ft)  for  this 
family  of  wave  spectra. 

It  is  noted  in  Figure  7  that  the  probability  of  breaking  waves  varies 
considerably  in  a  given  sea  severity  depending  on  the  shape  of  the  wave 
spectrum.  For  example,  in  a  sea  of  significant  wave  height  10.7  m  (35  ft), 
the  probability  ranges  from  3.97  percent  to  22.92  percent  with  a  most  probable 
of  10.38  percent  which  is  computed  by  using  the  most  probable  wave  spectrum 
of  the  family.  Figure  8  shows  three  spectra  of  the  family  for  the  sea  of 

significant  wave  height  10.7  m  (35  ft)  which  yield  the  largest,  most 


probable,  and  smallest  probability  of  breaking  identified  as  Spectra  I, 

II,  and  III,  respectively. 

The  results  of  computations  in  various  sea  severities  show  that  the 
frequency  of  occurrence  of  breaking  waves  is  a  function  of  the  4th  moment, 
m^  ,  of  the  wave  spectrum.  In  order  to  elaborate  on  this  statement.  Figure  (9) 
shows  the  probabilities  of  occurrence  of  breaking  in  six  different  sea 
states  of  the  six-parameter  wave  spectra  family  plotted  against  the  dimension¬ 
less  4th  moment  of  each  spectrum.  As  can  be  seen  in  the  figure,  the  frequency 
of  breaking  Increases  significantly  with  Increase  in  the  4th  moment  of  the 
wave  spectrum. 

The  conclusion  given  in  the  above  paragraph  can  be  verified  from  the 
breaking  condition  given  in  Equation  (4),  That  is,  from  the  definition  of 
the  average  time  interval,  T^  ,  Equation  (4)  can  be  expressed  in  the  following 
form: 


V  >  (4tt)^ 


-  (471)^  a 


l  +  (/l-e^)^ 


2^ 


(13) 


It  is  noted  that  /L  -e^/  ^1  +  A  -  j  in  the  above  equation  is  almost 
constant  for  the  practical  range  of  c-value  of  ocean  waves,  say  0.3  -  0.8. 
Therefore,  the  v-value  reduces  in  proportion  to  the  square  of  the  increase 
in  the  4th  moment,  m^  ,  of  the  wave  spectrum.  This  Implies  that  a  significant 
Increase  of  the  probability  of  breaking  is  expected  with  increase  in  the  4th 
moment. 
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In  order  to  substantiate  the  conclusion  given  above,  computations  of  the 
probability  of  occurrence  of  breaking  as  well  as  the  joint  probability  function 
of  wave  excursion  and  time  Interval  are  carried  out  for  four  spectra  measured 
in  the  North  Atlantic.  Two  of  the  spectra,  NW  1  and  NW  154  shown  in  Figure  10, 
are  those  in  a  sea  of  significant  wave  height  4.6  meters  (15.1  ft),  while 
the  other  two  spectra,  NU  314  and  JHC  75  shown  in  Figure  11,  are  those  in  a 
sea  of  10.8  meters  (35.4  ft).  The  results  of  the  computation  show  that 
there  is  almost  no  wave  breaking  for  the  two  spectra  in  the  sea  of  significant 
wave  height  4.6  meters  (1.1  percent  for  NW  1  and  none  for  NW  154).  On  the  other  hand, 
the  probability  of  breaking  is  11.90  percent  and  13.05  percent  for  the  spectrum  NW  134 
and  JHC  75,  respectively,  in  the  sea  of  significant  wave  height  10.8  meters.  A 
comparison  of  the  joint  probability  function  of  Type  I  excursion  and  time  interval, 
f(v>^)>  for  the  spectrum  NW  154  (significant  wave  height  4.6  m)  and  JHC  75  (sig-> 
nif leant  wave  height  10.8  m)  is  shown  in  Figure  12.  Included  also  in  the  figure 
is  the  line  indicating  the  breaking  condition.  As  can  be  seen  in  the  example 
shown  in  Figure  12,  the  shapes  of  the  dimensionless  joint  probability  density 
function  of  excursion  and  time  Interval  are  nearly  equal  for  these  sea  states. 

However,  the  location  of  the  line  indicating  the  breaking  condition  differs 
substantially  as  shown  in  the  figure  depending  on  the  magnitude  of  the  4th 
moment  of  the  spectrum.  This  in  turn  yields  the  difference  in  the  probability 
of  breaking  according  to  the  formulation  given  in  Equation  (13). 

It  is  of  Interest  to  note  that  the  relationship  between  the  probability 
of  breaking  and  the  4th  moment  of  the  spectrum  evaluated  for  various  spectra 
measured  in  the  North  Atlantic  agrees  very  well  with  the  functional  relation¬ 
ship  given  in  Figure  9  which  is  obtained  using  the  six-parameter  family  of 
wave  spectra. 

Since  breaking  increases  significantly  with  increase  in  the  4th  moment 
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of  the  wave  spectrum,  in  general,  breaking  is  very  sensitive  to  the  existence 
of  high  frequency  energy.  In  other  words,  breaking  is  more  likely  to  occur 
in  high  frequency  waves,  and  this  agrees  with  the  concept  of  the  equilibrium 
range  of  a  spectrum  (Phillips  1958). 

5.  Experiments  on  the  Frequency  of  Breaking  Waves 

An  experimental  study  was  carried  out  in  which  irregular  waves  were 
generated  in  the  tank  to  examine  whether  or  not  the  observed  frequencies  of 
occurrence  of  breaking  waves  agree  with  theoretical  results.  Unfortunately, 
due  to  the  limitation  in  the  capability  of  the  wave-maker,  it  was  difficult 
to  generate  high  frequency-component  waves  which  were  sufficiently  steep 
to  produce  breaking  associated  with  Type  n  excursions.  Although  records 
taken  during  the  tests  showed  a  number  of  Type  n  excursions  in  the  generated 
waves,  their  steepness  was  insufficient  to  result  in  breaking.  This  is  in  agree¬ 
ment  with  the  results  of  computations  carried  out  on  the  generated  waves  which 
indicate  the  probability  of  breaking  of  Type  II  excursions  is  extremely  small. 

In  order  to  compare  the  computed  and  observed  frequencies  of  wave 
breaking,  two  wave  probes  were  Installed  in  the  tank  separated  by  a 
distance  of  14.6  m  (48  ft).  An  example  of  wave  spectra  measured  at  these 
two  locations  is  shown  in  Figure  13.  As  can  be  seen  in  the  figure,  wave 
energy  was  substantially  reduced  due  to  wave  breaking  which  took  place 
while  the  waves  were  travelling  between  the  two  locations,  and  this  results 
in  a  reduction  in  the  magnitude  of  significant  wave  height  from  8.96  cm  to 
7.32  cm.  The  probability  of  breaking  waves  at  each  location  was  computed 
by  using  the  spectrum  measured  at  that  location  The  probability  of  breaking 
waves  thusly  computed  was  compared  with  that  obtained  experimentally.  The 
latter  was  evaluated  by  dividing  the  number  of  waves  visually  observed  to 
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break  between  these  two  locations  during  the  tests  by  the  total  number  of 
waves.  The  results  of  the  comparisons  made  for  four  (4)  different  sea 
conditions  generated  In  the  tank  are  tabulated  below  : 


Sea  condition 


n 

III 

IV 

0.015 

0.051 

O.ObO 

0.021 

0.043 

0.045 

Experimentally  obtained 
probability  of  breaking 


Computed  probability  of 
breaking 


As  can  be  seen  In  the  above  table,  reasonably  good  agreement  Is  obtained 
between  the  theoretical  and  experimentally  obtained  results. 

6.  Conclusions 

This  paper  discusses  a  method  to  evaluate  the  frequency  of  occurrence 
of  breaking  waves  In  deep  water.  Breaking  waves  defined  In  the  present 
study  are  those  associated  with  steepness,  and  the  breaking  criterion  was 
obtained  from  measurements  of  Irregular  waves  generated  In  the  tank.  It 
was  found  that  the  functional  relationship  between  wave  height  and  period 
at  the  time  of  breaking  of  the  Irregular  waves  Is  by  and  large  different 
from  that  known  for  regular  waves.  From  the  results  of  the  experimental 
study,  the  following  relationship  Is  obtained  as  the  criterion  for  breaking 
of  Irregular  waves: 

H  >  0.020  gT^ 

In  the  derivation  of  the  prediction  formula,  the  joint  probability 
distribution  of  wave  excursion  and  associated  time  Interval  for  a  non- 
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narrow-band  random  process  is  used.  Wave  breaking  which  takes  place  along 
an  excursion  crossing  the  zero-line  (Type  I  excursion)  as  well  as  that  which 
occurs  along  an  excursion  above  the  zero-line  (Type  n  excursion)  are  con¬ 
sidered. 

Comparisons  between  computed  and  observed  frequencies  of  occurrence  of 
wave  breaking  made  for  four  different  sea  conditions  generated  in  the 
tank  show  reasonably  good  agreement. 

Computations  of  the  probability  of  occurrence  of  breaking  are  made  in 
various  sea  severities  using  a  family  of  wave  spectra  developed  from  analysis 
of  wave  data  measured  in  the  North  Atlantic  Ocean.  From  results  of  compu¬ 
tations,  the  following  conclusions  may  be  drawn: 

(1)  In  evaluating  the  probability  of  breaking  waves,  the  frequencies  of 
occurrence  of  Type  I  and  Type  n  excursions  have  to  be  taken  into  consider¬ 
ation.  The  probability  of  breaking  waves  associated  with  Type  II  excursion 
cannot  be  neglected;  for  some  wave  spectra,  the  probability  is  approximately 
27  percent  of  that  associated  with  the  Type  I  excursion. 

(2)  The  results  of  computations  indicate  that  no  wave  breaking  is  expected  in  seas 

of  significant  wave  height  less  than  4  meters  (13  ft)  for  this  family  of  spectra. 

(3)  The  probability  of  occurrence  of  breaking  waves  depends  to  a  great 
extent  on  the  shape  of  the  wave  spectrum.  The  4th  moment  of  the  spectrisn 
is  a  dominant  parameter  which  Influences  the  occurrence  of  breaking.  The 
probability  Increases  significantly  with  increase  in  the  4th  moment  of  the 
spectrum  irrespective  of  sea  severity. 
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Appendix.  Joint  Probability  Density  Function  of  Type  n 
Excursion  and  Associated  Time  Interval 


In  order  to  derive  the  joint  probability  density  function  of  Type  H 
excursion  (BA  in  figure  1)  and  the  associated  time  interval,  T,  we  may  first 
consider  the  joint  probability  density  function  of  the  positive  minima,  B, 
and  the  associated  acceleration.  The  joint  probability  density  function 
will  then  be  transformed  to  the  joint  probability  density  function  of  the 
positive  minima  and  the  time  interval,  T. 

Let  N  be  the  expected  number  of  positive  minima  per  unit  time,  and 

D 

let  N-  ^  be  the  expected  number  of  positive  minima  per  unit  time  above  a 

^B,x  ^  ^ 

level  Cd  with  its  associated  acceleration  x.  These  are  given  as  follows 

D 

respectively , 


N..  = 


N, 


£  II 

^B,X 


^  jx  f(x,  o,  x)  dx  dx 
o  o 
00 

Jl»  II. 

X  r(x,  o,  x)  dx 


(A.l) 


f  If 

where,  f(x,x,x)  =  the  joint  probability  density  function  of  the 

displacement,  velocity  and  acceleration  of  a 
Gaussian  random  process  with  zero  mean  and  the 


covariance  matrix. 


Then,  the  joint  probability  density  function  of  the  positive  minima,  ? 


and  the  associated  acceleration,  x,  is  given  by. 


N-  •• 
1 - ^ 


N 


B 


X  f (Cg,  o,  x) 


(A.  2) 


w  oo 

If/ 


It  It 

f (x,  o,  x)  dx  dx 


*0  o 


By  carrying  out  the  integration  involved  in  the  denominator  of  (A. 2), 
we  have. 


fCCg.  X)  = 


II 

X 


^  (l-/l-  e^) 


(A. 3) 


0<Cg<“,  0<x<'» 


where,  A  =  m  m  -  m^ 

0  >•  2 

Next,  consider  the  joint  probability  density  function  of  the  positive 

ft 

maxima,  ^  ,  the  positive  minima,  ^  ,  and  associated  acceleration,  x. 


Since  it  may  be  safely  assumed  that  the  positive  maxima  and  the  positive 
minima  are  statistically  independent,  the  joint  probability  density  function 
is  given  simply  by  the  product  of  two  probability  density  functions, 

II 

f(  ^  )  and  f (  Cb  >  x) .  Here,  the  probability  density  function  of  the  positive 
A  D 

maxima  is  given  by  [Ochi  1973], 


Hence,  from  (A. 3)  and  (A. 4),  we  have. 


In  order  to  facilitate  the  integration  involved  in  further  analysis, 
the  following  approximation  is  made  in  equation  (A. 5); 
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Next,  let  us  consider  the  excursion  BA,  denoted  by  C  =  and 

derive  the  joint  probability  density  function  of  ,  C  »  and  x.  This 
can  be  simply  obtained  by  substituting  (C^-C)  for  in  equation  (A. 8). 
Then,  the  joint  probability  density  function  of  the  excursion  C  and  the 
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acceleration,  x,  can  be  obtained  as  the  marginal  probability  density  function 
of  f(C^  ,C  >  x).  That  is, 

ioo  OC 

•  fa^-^ ,  x)dc^  +  J  f(c^)  .  f(c^-c,x)dc^  (a.9) 

^  ^  <  oOj  0  <  X  <  «> 

A  U  ■* 


The  first  term  of  equation  (A.9)  is  for  C  <  0,  and  the  second  term  is 
for  ^  >  0.  Since  wave  breaking  is  associated  with  C.  >  t  namely  C  >  0, 

A  i) 

we  consider  only  the  second  term  for  the  present  problem.  The  integration 
of  the  second  term  of  equation  (A.9)  yields, 

j  f *  f  ,  x)  dC^ 


0  <  ^  <  00  0<x<“ 


It  is  noted  that  equation  (A. 10)  is  derived  by  truncating  the  probability 
density  function  given  in  (A. 9)  at  C  “  0*  Hence,  equation  (A. 10)  has  to  be 
multiplied  by  a  constant  such  that  it  satisfies  the  condition  required  to 
being  a  probability  density  function.  Taking  this  condition  into  consideration 
equation  (A. 10)  can  be  expressed  in  the  following  dimensionless  form: 


where,  v  ■  dimensionless  excursion  *  C/t^n~  with  C  =  -  Cb 

o  A  B 

M 

T  ■  dimensionless  acceleration  *  x//m 

>< 

The  function  k(e)  in  the  above  equation  is  associated  with  the  truncation 
of  the  probability  density  function  at  C  ■  0,  and  it  is  given  by. 


The  results  of  computations  using  equation  (A. 12)  for  various  c-values 
have  shown  that  k(e)  can  be  approximately  expressed  by  the  following  function 
of  E  with  sufficient  accuracy  in  the  practical  range  of  E-values  for  wind¬ 
generated  seas,  say,  0.3  <  E  <0.9: 


k(£)  =  1.13  +  0.011  e 


6.68(£  -  0.4) 


(A. 13) 


Next,  the  joint  probability  density  function  given  in  equation  (A. 11) 
can  be  converted  from  an  acceleration  to  a  time  interval  by  the  same  procedure 
as  was  used  for  the  derivation  of  the  joint  probability  density  function  of 
Type  I  excursion  and  the  associated  time  Interval  (Arhan,  M.  et  al.  1976, 
Cavanie,  A.  et  al.  1976,  Ezraty,  R,  et  al.  1977). 


It 

X 


(A. 14) 


In  the  dimensionless  form,  equation  (A. 14)  yields. 
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NOMENCLATURE 


i 


a 


f 

f(x,  X,  x) 


f(C.  X) 


f(CA. 


fCCg,  X) 


f(CA) 

fj(v,X),  fjj.(v,X) 


H 

L. 

m. 


Na.  Ng 


N 


BA 


wave  frequency  in  Hz 

joint  probability  function  of  the  displacement, 
velocity  and  acceleration  of  a  Gaussian  random  process 
joint  probability  function  of  the  excursion  and  the 
associated  acceleration 

joint  probability  density  function  of  the  positive 
maxima,  the  positive  minima  and  the  associated 
acceleration 

joint  probability  function  of  the  positive  minima  and 
the  associated  acceleration 

probability  density  function  of  the  positive  maxima 

joint  probability  density  function  for  the  dimensionless 
Type  I  or  Type  IE  excursion  and  the  associated  time  interval 
gravitational  acceleration 
wave  height 

Stokes  limiting  wave  length 

the  j-th  moment  of  the  wave  spectral  density  function 

expected  number  of  the  positive  maxima  and  minima  per 
unit  time,  respectively 

expected  number  of  Type  I  and  Type  IE  excursion  per  unit 
time,  respectively 


the  expected  number  of  positive  minima  per  unit  time 
above  a  specified  level  ?  with  its  associated  acceleration 

D 

•  I 

X 
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Figure  1  Explanatory  sketch  of  a  non-narrow-band 
random  process 


Figure  2  Examples  of  Irregular  waves  Imminent  to 


SPECTRAL  DENSITY  IN  M'^-SEC 


Figure  6  Wave  spectrum  used  for  computation  of  the  joint 
probability  density  functions  shown  in  Figures 
4  and  5.  Significant  wave  height  12.2  m  (40  ft) 
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Figure  10  Spectra  of  significant  wave  height  4,6 
measured  in  the  North  Atlantic 
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Figure  13  Examples  of  spectra  measured  at  two  locations 
In  the  tank 
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